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We demonstrate that 'Bimaximal + Democratic' type neutrino mass matrix can accommodate 
the deviation of 8q from its maximal value which referred in the literature as 'quark-lepton com- 
plementarity' along with the other present day neutrino experimental results, namely, atmospheric, 
CHOOZ, neutrinoless double beta decay (/3/3o v ) and result obtained from WMAP experiment. We 
define a function \ v in terms of solar and atmospheric neutrino mass squared differences and solar 
neutrino mixing angle (obtained from different experiments and our proposed texture) . The masses 
and mixing angles are expressed in terms of three parameters in our proposed texture. The allowed 
region of the texture parameters is obtained through minimization of the above function. The pro- 
posed texture crucially depends on the value of the experimental results of j3Po v experiment among 
all other above mentioned experiments. If, in future, Pf3o v experiments, namely, MOON, EXO, 
GENIUS shift the lower bound on (m ee ) at the higher side by one order, the present texture will be 
ruled out. 

PACS numbers: PACS number(s): 13.38.Dg, 13.35.-r, 14.60 .-z, 14.60.Pq. 



INTRODUCTION 



The conjecture of neutrino flavor oscillation has been 
strengthened by the Super-Kamiokande (SK) atmo- 
spheric neutrino experiment. This has also been substan- 
tiated by K2K long baseline neutrino experimental re- 
sults ■ The best-fitted values obtained from SK exper- 
iment are given by Am^ tm = 2.4 x 10~ 3 eV 2 , sin 2 29 atm 
= 1.0 3]. Moreover, global analysis of all solar neutrino 
experimental results , including SNO salt phase experi- 
ment, j3| is in favor of Large Mixing Angle (LMA) MSW 
solution of solar v e deficit problem and this has also in 
concordance with the KamLAND [5j reactor neutrino ex- 
perimental results. Furthermore, the CHOOZ Q exper- 
iment has also put an upper bound on the U e 3 element 
of the lepton mixing matrix and combined analysis of so- 
lar, atmospheric and reactor data put bound on U e s as 
|f/ e3 | < 0.22 at 99.73% c.l 0. There are two other experi- 
mental results, one of them, from WMAP experiment Q 
on cosmic microwave background anisotropics gives an 
upper bound on the total neutrino mass as Em, < 0.63 
eV 0. The second one, the neutrinoless double beta 
experiment |l0j| has put a bound on Major ana- type neu- 
trino mass as < m ee > = (0.05 - 0.84) eV at 99% c.l with 
an uncertanity of the nuclear matrix elements up to 50%. 

In view of the above it is worthwhile to investigate ap- 
propriate texture of lepton mass matrix which satisfies 
all those experimental results. First of all, exact 'bimax- 
imal' texture is not admissible due to departure of #0 
from 7r/4 and the deviation is parametrised in terms of 
Cabibbo angle and the phenomena is referred as 'quark- 
lepton complementarity' ^lj. Keeping the charged lep- 
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ton mass matrix flavor diagonal, in the present work, we 
consider a neutrino mass matrix of the form 'bimaximal 
+ democratic' where the 'democratic' part of the neu- 
trino mass matrix parametrise the difference 7r/4 — 9q. 
We have not addressed any typical model in this work, 
because, there are plenty of literature concerning bimaxi- 
mal structure of neutrino mass matrix , however, we have 
indicated possible source of democratic type matrix ex- 
plicit realization of which needs some extra flavor symme- 
try. The present texture gives rise to 9 a tm maximal and 
vanishing value #13. We fit two neutrino mass-squared 
differences with the required values of solar and atmo- 
spheric neutrino experimental results and the remaining 
mixing angle with the required values of solar neutrino 
neutrino mixing angle. These three parameters are fixed 
by defining a function \p (see later) as the sum of squares 
of the differences between the calculated values of neu- 
trino oscillation parameters (with the texture considered) 
and the best fitted values of the same (obtained from dif- 
ferent analyses) and then minimising this function. The 
plan of paper is as follows: In Section II, we propose the 
texture of the neutrino mass matrix and its eigenvalues 
and mixing angles. We discuss phenomenological fits of 
the model parameters through our defined \p function in 
Section III. Section IV contains our conclusion. 



II. PROPOSED TEXTURE 

We consider the following texture of the neutrino mass 
matrix 



M u = M u + M demo 
where the first part Mu is given by 
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and the second part Mdemo is as 

/I 1 1\ 

M demo = d 1 1 1 (2.3) 

V 1 1 V 

where d is some scale factor. Thus the entire mass ma- 
trix contains three parameters a, c and d and we consider 
all of them to be real. A possible source of the terms 
which gives rise to 'democratic' type mass matrix is due 
to higher dimensional operator 

i,j 

through considering equality between the /y couplings. 
Diagonalizing the entire neutrino mass matrix we obtain 
the following mixing angles 

023 = -t/4,031 = (2.5) 

tan 2 12 = L - 2.6 

m2 — a 

and the eigenvalues are 

3d + c + x 



3d + c — x 
m 2 = 

m 3 = -c (2.7) 

where 

x = ^(c + d) 2 + 8(a + df (2.8) 

where a = a' + d and c = d + d. Next, we set the solar 
and atmospheric neutrino mass squared differences as 

Am 2 oi = Amji = m\ — m 2 = x(3d + c) (2-9) 

and 

Am 2 tm = Am 2 23 = m 2 -m 2 = (3(d+c) -x){c-3d+x)/A 

(2.10) 

The best fit values of oscillation parameters from solar 
neutrino experiment , atmospheric neutrino experiments 
and solar neutrino mixing angle are used to obtain the 
values of a, c and d. Before going to the neumerical fit of 
the parameters, we would like to mention that the present 
mass matrix contains three parameters which are fitted 
with three experimental values, hence, the the texture 
has predictions about the neutrinoless double beta de- 
cay result and cosmological constraint on total neutrino 
mass (WMAP experimental result) . The predictibility of 
the texture would increase if we incorporate more sym- 
metries in an explicit model, which necessarilly reduce 



the number of parameters. However, since the present 
texture admits nicely the constraint from WMAP exper- 
iment ,the only testability of the present texture relies 
crucially on the future experimental results of neutrino- 
less double beta decay. If this experimental result gives 
a strong evidence in favour of Majorana-type neutrino 
mass, then the present texture should explain the result 
or it will be ruled out. In the next section we will fit 
those values with the best fit results through our defined 
Xp function. 

III. NEUTRINO PHENOMENOLOGY 

The best fit values of oscillation parameters from solar 
neutrino experiment and atmospheric neutrino experi- 
ments are used to obtain the values of a, c and d. For this 
purpose, we consider Am 2 2 , the difference of the square 
of mass eigenstates mi and m 2 and the mixing angle 
9i2 that are responsible for solar neutrino oscillation and 
Am 23 , 623 responsible for oscillation of atmospheric neu- 
trinos. A recent global analysis |4| of the solar neutrino 
data from all solar neutrino experiments namely Chlo- 
rine, Gallium, Super-Kamiokande, SNO charged current, 
neutral current including salt phase data, shows that 
large mixing angle or LMA solution is most favoured for 
solar neutrino oscillation. From the combined analysis 
with solar neutrino data and KamLand experiment the 
best fit values of the solar neutrino oscillation parameters 
are given by Am (= Am 2 2 for our model) = 8 x 10 -5 eV 2 
and tan 2 Q (= tan 2 $i 2 for our model) = 0.45 |4|. From 
the analysis of SK atmospheric neutrino oscillation data, 
we have the best fit values of Am 2 tm (= Am| 3 for our 
model) = 2.4 x 10~ 3 eV 2 and 6* atm is maximal. This value 
of # a tm has already been obtained in our proposed tex- 
ture for 823- Thus treating a, c and d as parameters we 
can obtain different values of Am 2 2 , Am 23 and tan 2 6*i 2 
and compare them with best fit values of those quan- 
tities namely Anig, Am 2 tm and tan 2 # obtained from 
the analysis of solar and atmospheric neutrino data (dis- 
cussed above) to fix a, c and d. To this end we define a 
function 

X 2 = (1 - Am 2 2 /Am|) 2 + (1 - Am 2 3 /Am 2 tm ) 2 

+ (1 - tan 2 0i 2 / tan 2 6» ) 2 . (2.11) 

The function x 2 as defined above is calculated for a wide 
range of values of a, c and d and the minimum of the 
function is obtained. The corresponding values of a, c 
and d are given below. 

a = -0. 145 eV 
c = -0. 152 eV 

d =0. 051 eV (2.12) 
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Am^, Am 2 3 and tan 2 6*12 obtained from the above values 
of a, c and d and their comparison with the best fit values 
for Atoq , Am 2 tm and tan 2 #© obtained from recent anal- 
ysis of the solar and atmospheric neutrino data are shown 
in Table 1. In order to find out the range of values of a, c 
and d that satisfy the la limits of Am 2 , and tan 2 O,-.-, for 
the (|4|) combined analysis of recent solar neutrino data 
and KamLand data, we have fixed the value of Am| 3 at 
its best fit value (shown in Table 1) and vary the pa- 
rameters a, c and d such that Am 2 2 and tan 2 812 satisfy 
the allowed range mentioned above. This range as given 
in Ref. is 7.6 x 10~ 5 eV 2 < Am| < 8.6 x lCr 5 eV 2 
and 0.38 < tan 2 8 Q < 0.54. The allowed region in pa- 
rameter space of a, c and d are shown in a 3D plot 
(Fig. 1). As seen from Fig. 1, the allowed region is 
confined in a very narrow band in the parameter space. 
We have also found the parameter space in a, c and 
d for which Amjg lies in the 90% c.l. range given by 
2.0 x 10- 3 eV 2 < Am|, < 5.0 x lO^eV 2 while the other 
oscillation parameters kept in the previous range. This 
range of Am| 3 with 90 % c.l for SK atmospheric neutrino 
data is obtained from Q . The allowed region in the pa- 
rameter space of a, c and d are shown in the 3D plot of 
Fig. 2. The value of d gives rise to the value of M as 

M = (<p) 2 /d ~ 10 15 GeV. (2.13) 

assuming (</>} = 250 GeV. This value of M could be 
recognized as the right-handed symmetry breaking scale 
in some models which contains an extra U(1)r or 
SU(2)r symmetry in which neutrino masses are gener- 
ated through see-saw mechanism. Again, the sum of the 
three neutrino mass comes out as Em, = 0.44 eV which 
is far below from the upper bound of WMAP experimen- 
tal result. Furthermore, the value of d obtained from the 



best fit result is at the lower end of the present limit on 
Majorana neutrino mass ((m ee ) = 0.05 — 0.84)eV 
obtained from the neutrinoless double beta decay experi- 
ment. Although the result is somewhat controversial [l^ . 
however, the testability of the present model lies cru- 
cially on the future experiments such as, EXO, MOON, 
GENIUS etc. 



IV. CONCLUSION 

We demonstrate that quark-lepton complementarity 
can be realised in a 'Bimaximal + democratic' type neu- 
trino mass matrix which includes two possible sources 
of neutrino mass. One of them is the usual 'Bimaxi- 
mal' type arise in different models and a possible source 
of 'Democratic' type is through the incorporation of 
higher dimensional mass terms assuming equlity between 
yukawa couplings. We then fix the parameters through 
our defined \p function with the two mass squared differ- 
ences solar and atmospheric and with the solar neutrino 
mixing angle. The texture admits a best fit result at the 
minimum value of x 2 . Furthermore, the texture gives 
rise to a slightly lower value of m ee mass compare to 
the experimental result of neutrinoless double beta de- 
cay experiment, and if the result survives in future with 
MOON, EXO, GENIUS experimental results, the present 
texture will be ruled out. 
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Present Work 


Experiment 


Amir, 
(eV 2 ) 


Ami, Rpf M 
(eV 2 ) 


o 09 v in -5 


cn y 1 n~ 5 

O.U A 1U 


(eV 2 ) 


AmL Ref.fJ 

(eV 2 ) 


2.6 x 1(T 3 


2.4 x 10~ 3 


tan^ 6*i2 


tan 2 Ref.[4] 


0.476 


0.45 



TABLE I: Best fitted values for neutrino oscillation parame- 
ters obtained in the present work. 



Figure Caption 



Fig. 1 The 3D plot showing the region of the parameters 
a, c and d that produce the values of At7i 2 2 and tan 2 #12 
within lex range of the best fitted values from global solar 
neutrino + KamLand data analysis Am| 3 remains 
fixed at the best fit value. See text for details. 



Fig. 2 Same as Fig. 1 but in this case instead of keeping 
Am^ fixed at the best fit value, it is varied within the 
90% C.L. range 2 x 1CT 3 < Am| 3 < 5 x lQ- 3 eV 2 0. 



0.054 
0.0535 

0.053 
0.0525 

0.052 
0.0515 

0.051 
0.0505 




-0.152 
0.153 
-0.154 
0.155 



d 



0.0545 r 

0.054 - 

0.0535 - 

0.053 - 

0.0525 - 

0.052 - 

0.0515 - 

0.051 - 

0.0505 - 



+ 



+ 



+ 



t 



+ + 



+ 



+ 



+ 



+ 
+ 



+ 



+ 



JP> -0.152 
^^0.154 
-0.156 



